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Abstract 

Important  ingredients  in  both  exact  and  approximative  treatments  that  aim  to  establish 
the  electromagnetic  response  of  an  individual  scatterer  are  so-called  depolarization  dyadics. 
Here,  the  focus  shall  be  on  both  the  derivation  of  closed-form  expressions  for  as  well  as 
numerical  evaluations  of  the  depolarization  dyadics.  It  will  be  shown  how  — from  the  initial 
approaches  to  the  subject  for  isotropic  mediums  — important  results  have  been  derived  for 
the  general  anisotropic  and  general  bianisotropic  regime  in  recent  years. 


1.  Introduction 

In  order  to  motivate  our  approach  to  the  topic  of  depolarization  dyadics,  we  shall  consider 
an  applicational  area  of  electromagnetics  in  which  they  play  a significant  role.  That  topic 
is  homogenization , where  two  or  more  component  mediums  (often  envisioned  in  particulate 
form)  are  mixed  together  to  form  a composite  material.  The  aim  of  homogenization  theories  is 
then  to  extract  the  constitutive  parameters  of  the  composite  material  from  a knowledge  of  the 
constitutive  and  geometrical  properties  of  the  constituent  mediums.  The  research  literature  on 
homogenization  is  vast.  The  interested  reader  will  find  many  important  historical  works  in  an 
anthology  [1]  and  a comprehensive  review  of  new  developments  in  the  homogenization  of  linear 
bianisotropic  materials  in  a recent  book  chapter  [2]. 

The  first  step  towards  the  derivation/estimation  of  the  composite’s  constitutive  properties 
requires  the  solution  of  an  electromagnetic  scattering  problem  in  which  the  electromagnetic  re- 
sponse of  an  individual  scatterer  must  be  determined.  This  task  can  be  achieved  by  adopting 
direct  or  indirect  scattering  approaches.  For  a detailed  discussion  of  merits,  difficulties  and 
limitations  of  these  approaches,  the  reader  is  referred  to  an  up-to-date  conceptual  review  by 
Lakhtakia  [3].  While  the  direct  scattering  approach  usually  consists  of  numerical  implementa- 
tions, some  analytical  progress  can  often  be  made  in  the  indirect  scattering  approaches,  especially 
when  certain  approximations  are  introduced  where  appropriate.  In  particular,  in  the  so-called 
Rayleigh  estimate , the  complete  information  about  the  electromagnetic  response  is  contained  in 
the  depolarization  dyadics. 

In  the  present  review  the  topic  of  depolarization  dyadics  is  developed  through  a complete 
and  exact  representation  of  the  electromagnetic  field  that  is  applicable  even  when  the  field  point 
is  contained  within  the  source  region.  Alternatively,  it  is  shown  how  the  depolarization  dyadics 
of  a general,  linear  bianisotropic  medium  can  be  extracted  with  a Fourier  technique  in  the  form 
of  surface  integrals.  Also,  a detailed  discussion  focusses  on  the  types  of  mediums  for  which  these 
integrals  can  be  explicitly  evaluated.  The  basis  of  this  review  on  depolarization  dyadics  is  a 
recent  review  paper  on  fields  in  the  source  region  [4]  which  may  be  consulted  for  further  details. 
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2.  Theoretical  Background 

Time-harmonic  electromagnetic  fields  are  being  considered  and  a time-dependence  of  exp(— icot) 
of  field  quantities  is  assumed  throughout  and  suppressed  (w  is  the  angular  frequency).  In  6- 
vector  notation,  the  complex-valued  field  and  source  phasors  are  defined  by1 


E(s)  = 


( Mix) 
{H(x) 


(1) 


where  E(x)  and  H{x)  are  standard  electric  and  magnetic  field  phasors  and  Je{x)  and  Jm(r) 
are  the  electric  and  magnetic  current  densities,  respectively.  The  most  general  linear  medium, 
commonly  referred  to  as  a bianisotropic  medium , is  described  by  the  frequency-dependent  con- 
stitutive relations 

(fg)=B-Efe>>  K=(!!).  <2> 

where  D{x ) is  the  dielectric  displacement  and  B(x)  is  the  magnetic  induction.  As  for  the  3x3 
constituent  dyadics  of  , g and  p are  the  well-known  permittivity  and  permeability  dyadics, 
whereas  £ and  £ fire  the  two  magnetoelectric  dyadics.  The  medium  parameters  will  normally  be 
complex  scalars  or  pseudoscalars  and  it  will  be  assumed  that  they  fulfil  the  general  covariance 
constraint:  Trace  ^£  • p~l  + pTl  * £)  = 0,  [5],  dictated  by  the  structure  of  modem  electromag- 
netic theory. 

Consequently,  the  Maxwell  equations  can  be  given  in  the  compact  form: 


[L(V)  + ia;K]  * E(x)  = Q(x) , 


L(V)  = 


a vx|\ 

-V  x I Q j ’ 


(3) 


where  f_  is  the  3x3  unit  dyadic.  Because  the  differential  operator  L(V)  and  the  constitutive 
dyadic  K are  both  linear,  (3)  permits  a solution  representation  of  the  form 

F(r)  = Fh{x)  + f G(x,  x')  • Q{x')d3x' , (4) 

Jv'  — — 

wherein  V'  is  a volume  such  that  Q(x)  = 0 for  x £ V'.  The  field  F^r)  is  a solution  of  (3)  for 
Q(r)  = 0,  i.e.,  mathematically  speaking  it  is  the  complementary  function.  The  entity  G(x,x') 
is  generally  called  the  dyadic  Green  function,  abbreviated  DGF  henceforth.  It  contains  the 
standard  3x3  DGFs  G , G , G and  G . It  now  follows  from  (3)  and  (4)  that  the  DGF 
fulfils  the  dyadic  differential  equation 

[L(V)  + iu;K]  -Q{x,x')  = d{x-x')l,  (5) 

where  S(x  — x')  is  the  Dirac  delta  function  and  I is  the  unit  dyadic. 

The  field  representation  (4)  holds  for  all  x.  While  it  can  be  applied  directly  in  the  case 
x V1,  special  care  is  required  when  xEV'  because  of  singularities  in  the  integrand.  That  case 
is  of  special  interest  in  this  context  as  it  leads  to  the  concept  of  the  depolarization  dyadics. 


3.  Field  Representation  with  the  Fikioris  Technique 

Two  broad  approaches  for  dealing  with  the  problem  of  correctly  evaluating  the  integral  in  (4) 
for  x€V'  have  been  developed.  The  first  of  these  was  pioneered  by  Van  Bladel  [6,  7],  and  it  was 

'3-vectors  (6-vectors)  axe  in  normal  (bold)  face  and  underlined  whereas  3x3  dyadics  (6x6  dyadics)  are  in 
normal  (bold)  face  and  underlined  twice. 


A 


V 
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later  systematically  generalized  by  Yaghjian  [8].  Their  work  and  that  of  others  during  the  1980s 
relates  to  isotropic  dielectric-magnetic  mediums  only,  i.e.,  for  | = el,  £ = Q,  C = 2j  /£  = PL- 

More  recent  work,  based  upon  the  principles  of  electromagnetostatics,  reported  some  progress 
in  the  application  of  that  technique  to  uniaxial  [9]  and  general,  homogeneous  bianisotropic 
mediums  [10].  The  essence  of  the  Van  Bladel/Yaghjian  technique  is  that  the  integral  in  (4)  is 
treated  by  excluding  a small  convex  region  (often  called  exclusion  region  or  exclusion  volume) 
surrounding;  the  source  point  of  the  integration  region.  The  integral  can  then  be  evaluated  and, 
eventually,  the  linear  dimensions  of  the  exclusion  region  are  shrunk  to  zero  — a procedure  that 
effectively  amounts  to  principal  value  integration. 

An  alternative  approach  was  developed  by  Fikioris  [11]  in  which  the  volume  of  the  exclusion 
region  is  not  required  to  be  infinitesimally  small  in  size.  As  a consequence,  estimates  for  the 
electromagnetic  field  in  the  source  region  thus  obtained  differ  from  those  arising  from  the  Van 
Bladel/Yaghjian  technique.  In  particular,  the  Fikioris  technique  gives  rise  to  a dependence  on 
the  electrical  size  of  the  exclusion  region  in  addition  to  that  on  its  shape;  see  [12]  for  further 
information.  Thus,  the  Fikioris  technique  amounts  to  the  regularization  of  a divergent  integral 

[13] .  It  is  also  worth  noting  that  the  Fikioris  technique  simplifies  to  the  Van  Bladel/Yaghjian 
technique  when  the  exclusion  volume  is  made  infinitesimally  small. 

While  these  developments  concerning  the  Fikioris  technique  also  related  to  isotropic  medi- 
ums only,  recent  years  have  seen  significant  extensions  of  the  method  to  more  complicated, 
linear  homogeneous  mediums.  Initial  application  of  this  approach  for  a bi-isotropic  medium 

[14]  was  followed  by  an  extension  to  a simple  uniaxial  bianisotropic  medium  (i.e.,  an  affinely 
transformed  bi-isotropic  medium)  [15].  The  Fikioris  technique  was  also  applied  to  an  uniax- 
ial dielectric-magnetic  medium  [16],  an  axially  uniaxial  bianisotropic  medium  (AUBM)  [17],  a 
simple  symmetric  bianisotropic  medium  [18]  and  an  affinely  transformable  AUBM  [19].  Results 
may  be  extended  to  mediums  related  to  these  by  affine  or  other  field  transformations  [20]. 

The  essence  of  the  Fikioris  technique  can  be  outlined  in  a straightforward  way.  In  order  to 
transform  (4)  into  a suitable  format,  a volume  V"  is  defined  in  such  a way  that  it  constitutes  a 
convex  region  surrounding  x — x'.  Also,  an  auxiliary  dyadic  Green  function  G_p(x,x'),  called 
the  Poisson  dyadic , must  be  introduced.  Hence,  (4)  may  be  rewritten  as 

F(x)  = Fh(x)  + f G(x,x')  • Q(x')  d3x'  + 

Jv'-v"  — — 

J [fife  *0  * Q(x')  — Q_P(x,  x')  • Q(x)  j (fix'  + g(x)  • Q(x) . (6) 

This  representation  is  an  exact  formula  without  any  simplifications,  applicable  to  a general, 
linear  bianisotropic  medium  as  described  by  the  constitutive  dyadic  in  (2).  In  connection  with 
the  field  representation  (6)  it  must  be  observed  that  an  explicit  formula  for  the  infinite-medium 
DGF  Q{x,x')  — and  consequently  also  for  the  Poisson  dyadic  Q.p(x,x')  — of  a general,  linear, 
homogeneous  bianisotropic  medium  is  not  available  at  this  time  (see  [21,  22]  for  review  works  on 
infinite-medium  DGFs  for  complex  mediums).  This  entails  then  that  the  dyadic  function  D(x) 
in  (6),  which  is  called  the  depolarization  dyadic , and  which  is  given  by 

D(x)  = [ G p(x,  x7)  (fix' , (7) 

— Jv"  — 

can  also  not  be  derived  for  a general,  linear  bianisotropic  medium. 

While  the  foregoing  developments  have  fully  defined  the  depolarization  dyadic  |2(x),  it  is  in 
approximations  of  (6)  that  its  role  and  significance  receives  further  illumination.  One  may  as- 
sume, for  example,  that  V'  is  itself  convex  in  shape  and  also  electrically  small,  i.e.,  its  maximum 
chord  length  is  smaller  than  a tenth  of  the  principal  wavelengths  in  the  medium,  say.  Then,  as 
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the  Fikioris  technique  does  not  require  V"  to  be  infinitesimally  small,  one  can  set  V'  = V".  As 
a consequence,  the  representation  of  the  field  simplifies  to  the  extent  that  the  integral  over  the 
region  V'  — V"  on  the  right-hand  side  of  (6)  vanishes. 

If  one  makes  the  additional  assumption  that  the  current  density  distribution  Q(x)  is  uniform 
inside  V'  = V",  one  obtains  the  long-wavelength  estimate  from  the  Fikioris  technique  as 


for  x € V',  where 


ILiwiz) -Eh fe)  + Mte)  + D(x)  • Q(x) , 


M(x)  = J [G_(x,  x')  — Q_p(x,  x7)j  (fix' . 


(8) 

(9) 


If,  finally,  the  term  involving  M(x)  is  ignored,  the  Rayleigh  estimate  of  the  electromagnetic  field 
in  the  source  region  arises  in  the  form 


Eflayfe)  - Efcfe)  + S(^)  * Q(s)  • (10) 

In  this  last  relation  we  can  now  directly  recognize  the  significance  of  the  depolarization  dyadic 
as  a mapping  from  the  source  Q(x)  to  the  scattered  field  F Boy(x)  — Fh(x).  Thus,  the  complete 
response  of  the  individual  scatterer  is  seen  to  be  contained  in  the  depolarization  dyadic  D(x). 


4.  The  Depolarization  Dyadic 

As  far  as  the  depolarization  dyadic  D(x)  is  concerned,  the  situation  appears  now  like  this:  if, 
for  a specific  type  of  medium  the  dyadic  Green  function  G(x,  x')  is  known,  one  can  obtain  the 
Poisson  dyadic  Q_p(x,x')  and  thus  finds  a representation  of  the  depolarization  dyadic  in  terms 
of  a volume  integral  by  virtue  of  (7).  Due  to  the  structure  of  G.p(x,  x'),  that  volume  integral 
can  always  be  transformed  into  a surface  integral  with  the  help  of  Gauss’  theorem. 

Consider,  for  example  (see  [16]  for  mathematical  details),  an  uniaxial,  dielectric-magnetic 
medium,  defined  by  constitutive  dyadics 


£ = e*I+  (ec  - et)  cc,  | = Q,  £=  Q,  = Vtl+ {v-c  ~ Vt) &Q, 

where  c is  a unit  vector.  After  convenient  decomposition  of  D(x)  according  to 

f D lx)  D (x)  \ 

D(x)  = [ =eeV-'  =smU  , 

= Id  (x)  d (x)  ) 

\ =mew  =mm v— y / 

into  3x3  component  dyadics,  we  find 

D =^-qeL  • e_1  , £>  =0,  D = 0,  D =—qmL  • uTl , 

— ee  ^ ^e  = e = 5 =em  =’  =me  =5  =mm  Hrn  =m  5 


where  qe  = ec/et  and  qm  — nd^t-  The  dyadics  m can  be  explicitly  evaluated  as 


=e,m  = \ ( ~~  + A(Qe,m)  ) (l  “ - A (qe,m)cc, 

AW  = T^(1-^7ptan','/?TT)' 


(11) 

(12) 

(13) 

(14) 

(15) 


A 


< 


V 


i 


However,  as  an  alternative  to  the  route  taken  thus  far,  it  turns  out  that  D(x)  may  be  obtained 
independently  through  a Fourier  space  analysis  which  does  not  require  a complete  knowledge  of 
the  infinite-medium  DGF  of  the  medium.  The  procedure  was  first  implemented  for  a general 
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anisotropic  dielectric  medium  [23]  (for  a spherical  exclusion  volume)  and  then  generalized  to 
the  most  general,  linear,  homogeneous  bianisotropic  medium  with  constitutive  dyadic  as  given 
by  (2)  (for  an  ellipsoidal  exclusion  volume)  [24].  Following  [24],  let  the  surface  of  an  ellipsoidal 
region  V'  be  parametrized  as  2^(0,  <f>)  — SR  • x(6,  <f>)  where  x(6,  (f>)  is  the  radial  unit  vector  in 
a spherical  coordinate  system  (with  angular  coordinates  6 and  (f>)  located  at  the  centre  of  V'\  S 
is  a linear  measure  of  the  size  of  V,  and  R is  the  symmetric  shape  dyadic  (R  = L when  V'  is 
spherical). 

Then,  the  3x3  depolarization  dyadics  in  (12)  are  given  by 

Rxy  = UTl  * Rxx,  • UTl , (A,  A'  = e,  m) . (16) 


The  3x3  dyadic  R^^  is  calculated  as  the  double  integral 


1 rn  (xmT...'x)xx 

= / # / dd  sin 0-jz — — — — — — — — — — , (17) 

_aa  Amuj  Jt=o  Je= o (®  • §!  • x)(x  • pf  • x)  — (x  • C * x)(x  • x)  v 

£ = url  • | • url , £'  = £_1  • i • url , c'  = url  • c • url , d = ur1  • a • url , (is) 


r = u , r — 

=ee  5 =em 


= -r 


T 


= -C'5  T =6'. 

^ 5 =mm  = 


(19) 


It  is  clear  that  in  general  the  double  integral  in  (17)  needs  to  be  evaluated  numerically. 
Explicit,  closed-form  evaluations  (for  ellipsoidal,  spherical,  cylindrical  and  cubical  exclusion 
volumes)  exist  for  the  isotropic  case:  both  for  achiral  isotropic  mediums  [8]  as  well  as  for  chiral 
mediums  [14].  For  uniaxial  dielectric-magnetic  and  for  uniaxial  bianisotropic  mediums  where 
at  least  one  of  the  medium  dyadics  g,  £,  C,  p is  of  the  form  g = acl  + accc  (c  is  again  a unit 
vector)  explicit  expressions  (in  terms  of  inverse  trigonometric  functions)  are  also  available  [23], 
[16],  [25];  [17],  [26].  The  most  recently  obtained  results  in  closed  form  relate  to  a biaxial  dielectric 
anisotropic  medium  with  g - exU^u^ . + eyUyUy  + ezuzuz  , £ = Q,  C = fl,  /£  = where  y^, 
■Uy,  Uj  are  the  unit  vectors  of  a cartesian  coordinate  system.  The  depolarization  dyadic  D(x) 
is  given  in  terms  of  elliptic  functions  of  the  first  and  second  kind  [27].  Finally,  the  closed-form 
results  obtained  for  isotropic/uniaxial/biaxial  mediums  can  be  trivially  extended  to  mediums 
which  contain  arbitrary  skew-symmetric  elements  in  their  3x3  constitutive  dyadics,  see  [24], 
[27],  because  any  skew-symmetry  is  filtered  out  in  the  integral  representation  (17). 


5.  Concluding  Remarks 

It  was  shown  that  depolarization  dyadics  play  an  important  role  in  electromagnetic  field  rep- 
resentations. They  arise  through  regularization  procedures  which  guarantee  that  the  field  rep- 
resentations also  remain  valid  when  the  field  point  is  in  the  source  region.  Their  significance 
can  be  interpreted  in  a straightforward  way,  namely  that  within  certain  well-defined  approxi- 
mations, particularly  the  Rayleigh  estimate,  they  provide  the  complete  information  about  the 
electromagnetic  response  of  an  individual  scatterer. 

The  depolarization  dyadics  of  a general,  linear  bianisotropic  medium  can  be  represented  in 
terms  of  a two-dimensional  integral.  There  are  many  special  cases  for  which  the  integrations  can 
be  performed  explicitly  to  yield  closed-form  expressions.  In  general,  however,  these  integrals  lend 
themselves  easily  to  numerical  treatment  that  require  only  small  computational  resources.  Such  a 
numerical  approach  is  also  economic  because  in  their  most  important  application,  depolarization 
dyadics  are  used  in  homogenization  calculations  that  are  necessarily  numerical  by  nature. 
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